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Abstract 

< 

y^ ' It is described how the Hermite-Pade polynomials corresponding to an algebraic 

approximant for a power series may be used to predict coefficients of the power series 

that have not been used to compute the Hermite-Pade polynomials. A recursive 

algorithm is derived and some numerical examples are given. 
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1 Introduction 

>: 

X 

^ ■ Using sequence transformation and extrapolation algorithms for the prediction 

of further sequence elements from a finite number of known sequence elements 
is a topic of growing importance in applied mathematics. For a short intro- 
duction see the book of Brezinski and Redivo Zaglia [3, Sec. 6.8]. We mention 
theoretical work on prediction properties of Pade approximants and related 
algorithms like the epsilon algorithm, and the iterated Aitken and Theta algo- 
rithms [2, 5, 7, 11], Levin-type sequence transformations [8, 9], the E algorithm 
[2, 10], and applications on perturbation series of physical problems [6, 8]. 

Here, we will concentrate on a different class of approximants, namely, the 
algebraic approximants. For a general introduction to these approximants and 
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the related Hermite-Pade polynomials see [1]. Programs for these approxi- 
mants are available [4] . We summarize those properties that are important for 
the following. 

Consider a function / of complex variable z with a known (formal) power 
series 

oo 



The Hermite-Pade polynomials (HPPs) corresponding to a certain algebraic 
approximant are A^ + 1 polynomials Pn{z) with degree p„ = deg{Pn), n = 0..N 
such that the order condition 

N 

Y^P^{z)f{zr = 0{z'^) (2) 

n=0 



holds for small z. Since one of the coefficients of the polynomials can be nor- 
malized to unity, the order condition (2) gives rise to a system of M linear 
equations for A^ + I]^=o Pn unknown polynomial coefficients. Thus, the coeffi- 
cient of z"^ of the Taylor expansion at z = of the left hand side of Eq. (2) 
must be zero for m = 0, . . . , M — 1. In order to have exactly as many equations 
as unknowns, we choose 

N 

M = N + J2pn (3) 

n=0 



and assume that the linear system (2) has a solution. Then, the HPPs Pn{z) 
are uniquely defined upon specifying the normalization. The algebraic approx- 
imant under consideration then is that pointwise solution a{z) of the algebraic 
equation 

N 

Po{z) + ^ Pn{z)a{zY = (4) 

n=l 



for which the Taylor series of a{z) coincides with the given power series at 
least up to order z^~^. 

We note that for A^ = 1, the algebraic approximants are nothing but the 
well-known Pade approximants. 

Although we assumed that the power series of / is known, quite often in 
practice, only a finite number of coefficients is really known. These coefficients 



then may be used to compute the Hermite-Pade polynomials and the algebraic 
approximant under consideration. 

We note that the higher coefficients of the Taylor series of a{z) may be con- 
sidered as predictions for the higher coefficients of the power series. The latter 
are also of interest in applications. 

The question then arises how to compute the Taylor series of a{z). If it is 
possible to solve the equation (4) explicitly, i.e. for A^ < 4, a computer algebra 
system may be used to do the job. But even then, a recursive algorithm for 
the computation of the coefficients of the Taylor series would be preferable in 
order to reduce computational efforts. 

In the following section, such a recursive algorithm is obtained. In a further 
section, we will present numerical examples. 



2 The recursive algorithm 

We consider the HPPs 

Pn 
Pfiiz) =J2Pn,jZ^ (5) 

i=o 
as known. Putting 

oo 

a{z) = J2akz'' (6) 

we obtain from Eq. (4) 

PO N Pn OO OO n 

E po.^' + E E Pn.z' E • • • E ^''-'■■■'-'" n «^^ = (7) 

j=0 n=lj=0 ki=0 fc„=0 m=l 



whence, by equating the coefficient of z'^ to zero, we obtain an infinite set of 
equations. Due to Eq. (2), all the equations for J < M are satisfied exactly 
fora, = /„j=0,...,M-l. 

As a ffist step, we compute om- We note that M > pq. Hence, the coefficient 
of z'^ does not involve any terms with poj. For this coefficient Rm, we only 
need to consider terms in Eq. (7) such that M = j + ki + . . . + kn and we 



obtain Rm = for 

N 



Rm — X! X! ^".i n ^kr, 

n=l j+ki+...+kn=M m=l 



The only terms on the RHS involving qm are obtained if exactly one of the 
km is equal to M, i.e., we have km = M, j = and kj = for j ^ m. Thus, 
we may rewrite all these terms as gmC where 

TV 

C=J2^Pn,0fr' (9) 

n=l 



and note that the rest Dm = Rm — (imC is independent of cm- Recalling 
Rm = 0, we obtain 

aM = -Dm/C (10) 



Proceeding analogously for J > M, only terms with J = j + ki + . . . + kn need 
to be considered. Hence, Rj = for 

N n 

Rj = Yl Y. Pn,j n ^k^ (11) 

n=l j+ki + ...+kn=J m=l 



Now, the only terms on the RHS involving aj are obtained if exactly one of 
the km is equal to J, i.e., we have km = J, j = and kj = for j ^ m. Thus, 
we may rewrite all these terms as ajC where C is defined above. Proceeding 
as before, we put Dj = Rj — ajC and obtain 

aj = -Dj/C (12) 



An equivalent form of the recursive algorithm is obtained in the following way: 
Consider for known P„ and ao, . . . , clj-i the expression 



U, 



J\dzJ 



N J 

E^n(^) E«.^'| (13) 

2=0 n=l \ jr=0 



It is easy to see, that this expression is exactly equal to Rj, and hence, is linear 
in the unknown aj. Thus, we may compute the quantities Dj by substituting 



aj = into Uj, which entails 



JldzJ 



N fj-l 

y: Pn{z) y: «. ^' I (14) 

z=0 n=l \j=0 



Thus, starting from J = M, one may compute all the aj consecutively by 
repeated use of Eqs. (9), (14), and (12). 

This concludes the derivation of the recursive algorithm. 



3 Modes of application 



Basically, there are two modes of application. 

a) one computes a sequence of HPPs and for the resulting algebraic approxi- 
mants, one predicts a fixed number of sofar unused coefficients, e.g., only one 
new coefficient. This mode is mainly for tests. 

b) one computes from all available coefficients certain HPPs. For the best 
HPPs one computes a larger number of predictions for sofar unused coeffi- 
cients. 

In the following examples, we concentrate on mode b). Here, it is to be ex- 
pected that the computed values have the larger errors the higher coefficients 
are predicted. 



4 Examples 



The examples serve to introduce to the approach. All numerical calculations 
in this section were done using Maple (Digits=16). 



Example 1 



As a first example, we consider N = 2, po = pi = p2 = 1 and, hence, M = 5. 
Since iV = 2, we are dealing with a quadratic algebraic approximant. Then, 
the recursive algorithm is started by aj = fj, j = 0, . . . , 4. For 05, we obtain 

Pl.1/4 +P2,1 (2/0/4 + 2/1/3 + /I) +P2,0(2/l/4 + 2/2/3) .^,. 

Pi,o + 2p2,o/o 



and for J > 5, we obtain 

Pl,iaj-1 + P2,l J2k=0 Oj-fc-lOfe + P2,0 J2k=l O-kO-J-k ,-. „x 

«J = -1^ 7 16 



For 

/(^) = (2 -3^)1/2 + 1/(5-^) (17) 

the HPPs are determined to be 

Po(^) = 1.- 1.544503593423590 z (18) 

Pi{z) = .1947992842134984+. 06783822675080703 z (19) 

P2{z) = -.5044536972622500-. 01090573365920830 z (20) 

The results for the predicted coefficients given in Table 1. 

Table 1 

The case oi N = 2, pq = pi = p2 = 1 foi Eq. (17). Displayed are the coefficients of 
the Taylor series, the predicted coefficients and absolute and relative errors of the 
predicted coefficients. 

j fj aj \fj-aj\ rel. error (%) 



5 


-.294 


-.294 


.001 


.18 


6 


-.330 


-.332 


.001 


.38 


7 


-.389 


-.392 


.002 


.58 


8 


-.475 


-.478 


.004 


.76 


9 


-.593 


-.599 


.006 


.93 


10 


-.756 


-.765 


.008 


1.10 



Example 2 



As a second example, we consider again N = 2, pq = pi = p2 = 1, and M = 5, 
but now the function 

/(^) = 17(l-2^)-V=^ + ^/(2-^) (21) 



with the HPPs 



Pq(z) = -49.52369318166839- 6.946105600281359 ;2 (22) 

Pi(^) = 1. + 1.695055482965655 2 (23) 

P2iz) = .1125387307324166 -.2732915349762758 2 (24) 



The results for the predicted coefficients given in Table 2. 

Table 2 

The case oi N = 2, po = pi = p2 = 1 for Eq. (21). Displayed are the coefficients of 
the Taylor series, the predicted coeffizients and absolute and relative errors of the 
predicted coefficients. 

j fj aj \fj-aj\ rel. error 



5 67.938 68.212 .274 .40 

6 120.739 122.291 1.552 1.29 

7 218.459 224.194 5.735 2.62 

8 400.498 418.053 17.555 4.38 

9 741.657 790.063 48.406 6.53 
10 1384.425 1509.437 125.012 9.03 



Example 3 

As a final example, we consider the case N = pq = pi = p2 = 2, whence 
M = 8, and the function 

f{z) = exp{z) (2-3 z)-'^/^ + 1/(5 - z) (25) 

The corresponding HPPs are 

Pq(z) = 1. - 1.027576803009053 z + .02070967420422950 z'^ (26) 

Pi{z) = 2.617867885747464 -.6563757889994458 2 

-3.118191126500581^2 (27) 

P^^z) = -3.647182626894738 + 7.471780741166546 z 

- 3.356878399103086 z^ (28) 

The results for the predicted coefficients are displayed in Table 3. 



Table 3 

The case oi N = 2, po = pi = p2 = 2 for Eq. (25). Displayed are the coefficients of 
the Taylor series, the predicted coeffizients and absolute and relative errors of the 
predicted coefficients. 

j fj aj \fj-aj\ rel. error (%) 



8 


3.888956 


3.878509 


.010447 


.27 


9 


5.356681 


5.301047 


.055634 


1.04 


10 


7.451679 


7.275227 


.176452 


2.37 


11 


10.447061 


10.006950 


.440111 


4.21 


12 


14.739132 


13.781978 


.957155 


6.49 


13 


20.903268 


18.995972 


1.907297 


9.12 



5 Conclusions 



It is seen that even rather low-order algebraic approximants, or HPPs, respec- 
tively, can lead to quite accurate predictions of the unknown coefficients of 
the power series, especially for /m, and the next few coefficients. 
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